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ABSTRACT 

Wc present a detailed auto-consistent study of the nearest blue compact dwarf galaxy 
NGC 6789 by means of optical and UV archive photometry data and optical long- 
slit ISIS-WHT spectroscopy observations of the five brightest star-forming knots. The 
analysis of the spectra in all knots allowed the derivation of ionic chemical abun- 
dances of oxygen, nitrogen, sulphur, argon and neon using measures of both the 
high- and low-excitation electron temperatures, leading to the conclusion that NGC 
6789 is chemically homogeneous with low values of the abundance of oxygen in the 
range 12-|-log(0/H) = 7.80-7.93, but presenting at the same time higher values of the 
nitrogen-to-oxygen ratio than expected for its metal regime. 

We used archival Hubble Space Tefescope/Wide-Field Planetary Camera 2 
(HST/WFPC2) F555W and F814W observations of NGC 6789 to perform a pho- 
tometric study of the colour-magnitude diagram (CMD) of the resolved stellar pop- 
ulations and derive its star formation history (SFH), which is compatible with the 
presence of different young and old stellar populations whose metallicities do not nec- 
essarily increase with age. We fit the observed optical spectrum in all the five knots 
using the STARLIGHT code and a combination of single stellar populations following 
the SFH obtained from the CMD. We compare the resulting stellar masses and the 
relative fractions of the ionising populations with a non-constrained SFH case. The 
properties of the younger populations were obtained using CLOUDY photoionisation 
models, giving similar ages in all the knots in the range 3-6 Myr and the estimation 
of the dust absorption factor, which correlates with the observed GALEX FUV-NUV 
colour indices. The total photometric extinction and dust-absorption corrected Ha 
fluxes were finally used to derive the star formation rates. 

Key words: ISM: abundances - H ii regions - galaxies: dwarf - galaxies: individual: 
NGC 6789 ~ galaxies: starburst - galaxies: star formation - Hertzsprung-Russell and 
colour-magnitude diagrams 



1 INTRODUCTION 

Blue compact dwarf galaxies (BCDs) are known by their in- 
tense processes of star-formation. These make that their op- 
tical spectra were dominated by the blue light from massive 
young stars and the bright emission lines from the ionised 
gas, so they are also called H ll galaxies. 

Their high specific star-formation rates and their low 
gas-phase metallicities were taken as evidences that these 
galaxies were undergoing their first bursts of star formation 
(Sargent & Searle 1970). In fact, one of the most important 
aspects about BCDs is that they constitute an important 
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link to the high redshift universe and the early epoch of 
galaxy formation (Bergvall & Ostlin 2002). During the last 
decade, however, deep imaging of BCDs showed that in ad- 
dition to bright young stars from the present starburst most 
of them have an underlying older stellar population (Pa- 
paderos et al. 1996; Aloisi et al. 1999; Ostlin 2000). 

One of the challenges of the BCDs is to reconcile their 
low-observed metallicity with the relatively high star forma- 
tion rate (SFR) in these galaxies. Matteucci & Chiosi (1983) 
proposed three possible mechanisms to explain this fact; 1) 
variations in the initial mass function (IMF); 2) accretion 
of metal-poor gas; and 3) galactic winds powered by super- 
novae explosions. According to numerical models, it seems 
that the last mechanism is the most simple to reproduce the 
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observed properties (see Tolstoy et al. 2009 and references 
therein for a detailed explanation). 

The star formation history (SFH) and the metal con- 
tent of BCDs must be figured out in order to shed some light 
on these questions. In the case of SFH, although some ad- 
vances have been done by means of the fitting of the optical 
spectrum with stellar populations synthesis (Perez-Montero 
et al. 2010), one of the most reliable method consists of the 
analysis of colour-magnitude diagrams (CMDs). However, 
this is limited to the closest objects whose stellar popula- 
tion can be resolved. Another important contribution comes 
from the study of the ionic chemical abundances of elements 
with a different nucleosynthetic origin, such as oxygen and 
nitrogen (MoUa et al. 2006). These chemical abundances 
can be derived very accurately in the metal-poor gas phase 
of BCDs, where the cooling rate is not efficient, increas- 
ing the electron temperature of the ionised gas and enhanc- 
ing the emissivity of the coUisional lines necessary to derive 
the ionic abundances following the method based on the de- 
termination of the electron temperature (Perez-Montero & 
Diaz 2003; Hagele et al. 2006). 

NGC 6789 is, according to Drozdovsky et al. (2001), the 
closest BCD to our Galaxy, with a distance modulus of (m 
- M) = 27.80, or a distance of 3.6 Mpc. This value yields 
a linear scale of 17.5 pc arcsec"^. Drozdovsky & Tikhonov 
(2000) imaged NGC 6789 from the ground. They found that 
it belongs to the iE subtype (following Loose & Thuan 1986 
classification), exhibiting the morphology most characteris- 
tic of the vast majority of BCDs. NGC 6789 presents several 
Ha emission knots which show evidence of actual star for- 
mation activity. It also shows a high surface brightness in its 
central region and a small radial velocity (Vo = -141 km s~^, 
Karachentsev & Makarov 1998). NGC 6789's closeness, spa- 
tial isolation and morphology offer the prospect of studying 
the structure, SFH, and metallicity of different star forming 
knots within the same BCD. 

To investigate these issues, we obtained simultaneous 
blue and red long-slit observations with the ISIS double- 
arm spectrograph at the 4.2m William Herschel Telescope 
(WHT), of the brightest five knots of NGC 6789. We also 
used archival data from HST/WFPC2, GALEX and Ha nar- 
row filter, and B and R broad filters. These images were used 
to derive the SFH and to give an observational input to de- 
rive the properties of the ionising populations with the aid 
of photoionisation models. These were also used to provide 
more accurate chemical abundances in each one of the five 
studied knots of this galaxy and to better understand their 
differential chemical evolution. 

In the following section the long-slit WHT observations 
and reduction are described and the results from the analysis 
are presented. In section 3 the optical and UV photometry 
are described together with the resolved stellar photometry. 
Finally, we discuss all these results in Section 4 and conclu- 
sions are presented in Section 5. 



mounted on the 4.2 m William Herschel Telescope (WHT) of 
the Isaac Newton group (ING) at the Roque de los Mucha- 
chos Observatory on the Spanish island of La Palma. They 
were acquired on 2005 July 8 during one single night observ- 
ing run and under photometric conditions, with an average 
seeing of 0.7 arcsec. The EEV12 and Marconi2 detectors 
were attached to the blue and red arms of the spectrograph, 
respectively. The R600B grating was used in the blue cov- 
ering the wavelength range 3670 - 5070 A (centred at Ac — 
4370 A), giving a spectral dispersion of 0.45 A pixeP^. On 
the red arm, the R316R grating was mounted in two dif- 
ferent central wavelengths providing a spectral range from 
5500 to 7800 A (Ac = 6650 A) and from 7600 to 9900 A 
(Ac = 8750 A) with a spectral dispersion of 0.86 A pixel"^. 
To reduce the readout noise of our images, the observations 
were taken with the "SLOW" CCD speed. The pixel size 
for this set-up is 0.2 arcsec for both spectral ranges. The 
slit width was 1 arcsec, which, combined with the spectral 
dispersions, yields spectral resolutions of about 1.0 and 3.5 
A FWHM in the blue and red arms, respectively. The in- 
strumental configuration and other details on the exposures 
are given in the journal of observations in Table 1. 

Several bias and sky flat fleld frames were taken at the 
beginning and at the end of the night in both arms. In addi- 
tion, two lamp flat flelds and one calibration lamp exposure 
were taken for each telescope position. The calibration lamp 
used was -l-CuAr. 

The images were processed and analysed with IRAF^ 
routines in the usual manner. The procedure included cos- 
mic rays removal, bias subtraction, division by a normalised 
flat field and wavelength calibration. Typical wavelength fits 
to second to third order polynomials were performed us- 
ing around 40 lines in the blue and 20-25 lines in the red. 
These fits were done at 100 different locations along the 
slit in both arms (beam size of 10 pixels) obtaining rms 
residuals between ~0.1 and ~0.2pix. In the last step, the 
spectra were corrected for atmospheric extinction and flux 
calibrated. For both arms, BD-(-254655 standard star ob- 
servations were used, allowing a good spectrophotometric 
calibration with an estimated accuracy of about 5%. 

Unfortunately, the sky subtraction in the spectral range 
from 7600 to 9900 A left strong residuals from night-sky 
emission lines and telluric absorption, leaving this part of the 
spectra unusable and therefore, we were not able to measure 
with enough accuracy the [Slll] AA 9069,9532 A lines in any 
of the five knots. 

Figure 1 shows the Ha image and continuum contours 
in the R-band from Gil de Paz et al. (2003) in order to 
illustrate the position of the slits in relation to the posi- 
tion of the bursts of star formation and the host galaxy. We 
also show elliptical regions in Ha taken to measure the total 
Ha and GALEX fiux of the knots (see Section 3.1 below). 
The five knots are labelled with the first alphabet letters. 
Knots A and C were taken at PA = 138°, which in aver- 
age is close to the parallactic angle. In the case of knots B, 



2 LONG-SLIT SPECTROSCOPY 

2.1 Observations and reduction 

The long-slit spectrophotometric observations of NGC 6789 
were obtained using the ISIS double-beam spectrograph 



^ IRAF: the Image Reduction and Analysis Facility is distributed 
by the National Optical Astronomy Observatories, which is op- 
erated by the Association of Universities for Research in Astron- 
omy, In. (AURA) under cooperative agreement with the National 
Science Foundation (NSF). 
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Table 1. WHT instrumental configuration 



Slit position 


Spectral range 

(A) 


Disp. 

(Apx-i) 


FWHM 

(A) 


Spatial 
("px- 


res. 


Exposure Time 

s 


SI 


3670-5070 


0.45 


1.0 


0.2 




4 X 900 


SI 


5500-7800 


0.86 


3.5 


0.2 




2 X 900, 1 X 300 


SI 


7600-9900 


0.86 


3.5 


0.2 




2 X 900 


S2 


3670-5070 


0.45 


1.0 


0.2 




5 X 900 


S2 


5500-7800 


0.86 


3.5 


0.2 




2 X 900, 1 X 300 


S2 


7600-9900 


0.86 


3.5 


0.2 




2 X 900 



63° 




00" - 



44 



43' 



19'' 16"" 42' 

R.A. (2000.0) 



41' 



40' 



39' 



Figure 1. Ha image of NGC 6789 from Gil de Paz et al. (2003) with the identification of the observed knots (labelled from A to E). We 
show the regions for which the Ho flux was measured and the position of the two slits for the observations described in the text. The 
R-band contours show the position of the host galaxy. Finally, we also plot part of the rectangle encompassing the field-of-view of the 
PC chip of the WFPC2 data used to obtain the resolved stellar photometry described in Section 3.2. North is up and cast is towards the 
left-hand side. 



D, and E the slit was at PA = 76°, which is in average « 
40° off the parallactic angle. Therefore, the observations in 
this second position can be affected by a certain differential 
atmospheric refraction (DAR). Nevertheless, taking into ac- 
count the curves given by Filippenko (1982) for the mean air 
mass during the observations (~ 1.30), we calculated that 
the angular deviation between [On] and [Sll] is not larger 
than the slit width at this P.A. Besides, we checked that 
the emission line-ratios in the spatial position where both 
slit positions intersect do not vary in more than the quoted 



errors. Nevertheless, the line measurements errors of the sec- 
ond slit position (B, D and E knots) might be larger due to 
this effect. For the sake of comparison, we also plot in Fig- 
ure 1 part of the rectangle encompassing the field-of-view 
(FOV) of the PC chip of the WFPC2 data used to obtain 
the resolved stellar photometry (see Section 3.2) 
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Figure 2. Blue and red spectra for the knots A, B, C, D, and E of NGC 6789. 
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Figure 2. — continued 



2.2 Line intensities and reddening correction 

The optical calibrated spectra of the five observed knots with 
some relevant identified emission lines are shown in Figure 
2. The spectrum of each knot is split into two panels. 

Underlying stellar population in star forming galaxies 
have several effects in the measure of the emission lines pro- 
duced by the ionised gas. Balmer and Paschen emission lines 
are depressed by the presence of absorption wings of stellar 
origin and does not allow the measure of their fluxes with 
acceptable accuracy (Diaz 1988). All the properties derived 
from ratios that involve these lines, like reddening or ionic 
abundances, will be affected. 

We subtracted from the observed spectra the spec- 
tral energy distribution of the underlying stellar popula- 
tion found by the spectral synthesis code STARLIGHT" 
(Cid Fernandes et al. 2004, 2005; Mateus et al. 2006). 
STARLIGHT fits an observed continuum spectral energy 
distribution using a combination of the synthesis spectra of 
different single stellar populations (SSPs; also known as in- 
stantaneous burst) using a x^ minimisation procedure. We 
chose for our analysis the SSP spectra from Bruzual & Char- 
lot (2003), based on the STELIB hbrary of Le Borgne et al. 
(2003), Padova 1994 evolutionary tracks, and a Chabrier 
(2003) IMF between 0.1 and 100 Mq. The metallicity of 
the libraries was strictly constrained following the results 
derived from Section 4.1. Thus, we fixed the metallicity of 
the stellar populations to Z = 0.004 (« 3/10 Z©'') for the 
interval log(age) = [8.1,8.5] and Z = 0.001 (^ 1/12 Z©) 
for log(age) = {6.0, 6.4, 6.6, 6.8, 7.0, 7.1, 7.2, 7.9, 8.0, 8.1, 
8.2, 8.3, 9.1, 9.2, 9.3, 9.7, 9.8, 9.9, 10.0, 10.1} shown in the 
SFH solution (see Figure 5). As explained in section 4.1, al- 
though both metallicities appear in the youngest burst (< 
12 Myr) SFH solution, we assigned the lower metallicity to 
this event according to the abundances derived in section 
2.4. The STARLIGHT code solves simultaneously the ages 
and relative contributions of the different SSPs and the aver- 
age reddening. The reddening law from Cardelli et al. (1989) 
with Rv ~ 3.1 was used. Prior to the fitting procedure, the 
spectra were shifted to the rest frame and re-sampled to a 



^ Tlie STARLIGHT project is supported by the Brazilian agen- 
cies CNPq, CAPES and FAPESP and by the France-Brazil 
CAPES/Cofecub program. 
3 Zq = 0.0122 



wavelength interval of 1 A in the entire wavelength range by 
interpolation, as required by the program. Bad pixels and 
emission lines were excluded from the final fits. 

It should be noted that while emission lines can be 
masked out, this is not possible for the nebular continuum 
emission. We tested the contribution of this emission using 
Starburst99 libraries (SB99; Leitherer et al. 1999) in the 
same STARLIGHT models and we checked that the nebular 
continuum does not affect significantly neither the subtrac- 
tion of the underlying continuum nor the determination of 
the stellar mass, consistently with the weak gas emission in 
this galaxy. 

The Hi series (emitted as a consequence of recombina- 
tion) were used to determine the extinction, comparing the 
observed line ratios with the expected theoretical values. 
Case B (optically thick in all the Lyman lines) is the best 
simple approximation to describe the physical conditions in 
the ionisation of the gas. This method takes advantage of 
the fact that the ratio between the emissivities of two hy- 
drogen recombination lines, which depends on electron tem- 
perature and density, is almost constant. As an example, 
the ratio between the emissivity of Ha and H/3 is 2.86 for 



the case B with n. 



100 cm--^ and Te = 10000 K, and 



this value varies less than 10% in the range of interest of 
temperatures and densities for an Hll region. We used an 
iterative method to estimate them, taking as starting values 
those derived from the measured [Sll] AA 6717,6731 A and 
[Qui] AA 4363,4959,5007 A. A least square fit of the mea- 
sured decrements to the theoretical ones, computed based on 
the data by Storey & Hummer (1995), was performed that 
provides the reddening coefficient, c(H/3), and adopting the 
extinction law given by Cardelli et al. (1989) with Rv = 
3.1. Due to the large error introduced by the presence of the 
underlying population, only the strongest Balmer emission 
lines (Hof, H/3, H7, and H(5) were used. 

Line fiuxes for the most relevant emission lines were 
measured using the splot task in IRAF. Balmer lines were 
obtained from the STARLIGHT residual spectra, while the 
rest of the lines were measured on the original spectra. We 
checked in all spectra that the Starlight models properly 
reproduce the observed Balmer absorption profiles^. The 



From several fitting experiments using other libraries it is found 
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statistical errors associated with the observed emission line 
fluxes were calculated using the expression: 



Ul — (Jc 



N{ 1 + 



EW \ 
NAJ 



(1) 



where ai is the error in the observed line flux, (Jc represents 
the standard deviation in a box near the measured emission 
line and stands for the error in the continuum placement, 
N is the number of pixels used in the measure of the line 
flux, EW is the line equivalent width, and A is the wave- 
length dispersion in A per pixel (Gonzalez-Delgado et al. 
1994). This expression takes into account the error in the 
continuum and the photon count statistics of the emission 
line. 

Table 2 gives the equivalent widths and the emission- 
line fluxes relative to 1000-F(H/3), before and after reddening 
correction, in the optical spectra of the flve observed knots, 
together with the reddening constants and their errors, their 
corresponding A(V) value'' and the extinction-corrected H/3 
flux. We also provide the adopted reddening curve, /(A) nor- 
malised to H/3. The errors in the emission-line ratios were 
obtained by propagating in quadrature the observational er- 
rors in the emission-line fluxes and the reddening constant 
uncertainties. 



2.3 Electron densities and temperatures 

The electron density and temperatures of the ionised gas 
were derived from the emission line data using the same 
procedures as in Perez-Montero & Diaz (2003), based on the 
flve-level statistical equilibrium atom approximation in the 
IRAF task temden (De Robertis et al. 1987; Shaw & Dufour 
1995). See Hagele et al. (2008) and Appendix B of Garci'a- 
Benito (2009) for a description of the temperature and den- 
sities relations, respectively. We took as sources of error 
the uncertainties associated with the reddening-corrected 
emission-line fluxes and we propagated them through our 
calculations. 

Electron densities were derived from the [Sll] AA 
6717,6731 A and [On] AA 3726,3729 A line ratios, which are 
representative of the low-excitation zone of the ionised gas. 
In all cases they provide upper limits, which are remarkably 
similar for both ions. In the case of the oxygen lines, the 
spectral dispersion does not allow the total resolution of the 
lines, and they were deblended by a multi-Gaussian fit. The 
upper limits for the electron density are lower in all cases 
than the critical value for coUisional de-excitation. 

We computed three electron temperatures, T([On]), 
T([Olll]) and T([Slll]), for each of the flve knots. The auroral 
[Olll] 4363 was detected in all the five knots with sufficient 
signal-to-noise, and the [Olll] electron temperature was de- 
rived directly from the ratio (1(4959 A) + 1(5007 A))/ 1(4363 

A). 

The auroral [On] 7319 and [Oil] 7330 lines from the 2F 



that the resolution of the SSP spectra is a critical factor for 

Starlight in order to fit the absorption fines. 

5 c(H,9) = 0.4656*A(V), using Cardelli et al. (1989) and Ry = 

3.1. 



multiplet are emission-line doublets. Due to the spectral res- 
olution, we report here the sum of each doublet. Using the 
calculated [Olll] electron temperatures, we checked that the 
contribution by direct recombination is negligible for these 
lines (Liu et al. 2000). For knot A and D, the intensity of 
the [On] AA 7319,7330 A did not allow an accurate measure, 
therefore we derived their [On] temperature from T([Oni]) 
using the relation based on the photoionisation models de- 
scribed in Perez-Montero & Diaz (2003), which takes into 
account explicitly the dependence of T([Oll]) on the elec- 
tron density. 

Although we could measure the [Sin] A 6312 A in some 
of the knots, we were not able to get the [Sin] A 9069,9532 
A lines, so we estimated the [Sin] temperature from the 
empirical relation 

i([Siii]) = (1.19 ± 0.08)f([Oiii]) - (0.32 ± 0.10) (2) 

found by Hagele et al. (2006). 

The electron densities and temperatures derived in the 
flve knots of NGC 6789 are listed in Table 3 along with their 
corresponding errors. 



2.4 Ionic abundances 

We derived the ionic abundances of the different chemical 
species using the brightest available emission lines detected 
in the optical spectra using the task ionic of the STSDAS 
package in IRAF (see Hagele et al. 2008). The total abun- 
dances were derived by taking into account, when required, 
the unseen ionisation stages of each element, resorting to 
the ionisation correction factors (ICFs) for each species de- 
rived from the photoionisation models described in Section 
4.2 below. 



_.JCF(X')_ 



(3) 



2.4-1 Helium abundance 



Four of the strongest helium emission lines, Hel A A 4471, 
5876, 6678, and 7065 A were used to calculate He-|-/H-|-. 
The last two lines were only detected with enough signal-to- 
noise in knot C. Also in this knot. He n A 4686 was measured, 
allowing the calculation of twice ionised He. 

Helium lines arise mainly from pure recombination; 
however, they could have some contribution from coUisional 
excitation as well as being affected by self-absorption and, 
if present, by underlying stellar absorptions (see Olive & 
Skillman 2001, for a complete treatment), but we checked 
that these absorptions produce negligible variations in the 
measure of the lines in the STARLIGHT subtracted spectra. 

We took the electron temperature of [Olll] as represen- 
tative of the zone where the He emission arises and we used 
the equations given by Olive & Skillman (2001) to derive 
the He"'"/H'^, using the theoretical emissivities scaled to H/3 
from Smits (1996) and the expressions for the coUisional 
correction factors from Kingdon & Ferland (1995). To cal- 
culate the helium twice ionised, we used the equation found 
by Kunth & Sargent (1983). We did not take into account, 
however, the corrections for fluorescence since the involved 
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Table 2. Observed and reddening corrected relative line intensities [F(H/3)=I(H/3)=1000] with their corresponding errors for the five 
knots. The adopted reddening curve, /(A) (normalised to H/3), the equivalent width of the emission lines, the extinction-corrected H/3 
intensity, the reddening constant c(H/3) and the corresponding A(V) are also given. 







WCC 6789-A 




WGO 6789-B 


A (A) 


f(A) 


F{A) 


EW{A) 


1(A) 


F(A) 


EW(A) 


1(A) 


3726 [On] 


0.322 


896 ± 19 


-15.0 


979 ± 63 


1186 ± 17 


-34.9 


1277 ± 80 


3729 [On] 


0.322 


1241 ± 16 


-20.7 


1357 ± 84 


1744 ± 16 


-51.9 


1878 ± 116 


3835 H9 


0.299 


37 ± 10 


-0.6 


40 ± 11 


80 ±4 


-3.5 


85 ± 7 


3868 [Neiii] 


0.291 


302 ± 17 


-2.8 


327 ± 26 


153 ± 14 


-3.7 


164 ± 17 


3889 Hei+H8 


0.286 


153 ± 6 


-2.6 


166 ± 12 


146 ± 10 


-4.5 


156 ± 14 


3968 [Neiii]+H7 


0.266 


81 ± 8 


-1.3 


88 ± 10 








4102 ns 


0.229 


247 ± 5 


-3.9 


263 ± 16 


248 ± 16 


-6.7 


261 ± 22 


4340 H7 


0.157 


454 ± 11 


-7.2 


474 ± 28 


448 ± 7 


-12.4 


465 ± 26 


4363 [Oiii] 


0.149 


63 ± 3 


-0.7 


66 ± 5 


23 ±4 


-0.4 


23 ± 4 


4471 Hei 


0.115 


36 ± 4 


-0.4 


37 ± 5 


37 ± 6 


-0.7 


39 ± 7 


4861 H^ 


0.000 


1000 ± 10 


-17.3 


1000 ± 47 


1000 ± 12 


-29.1 


1000 ± 48 


4959 [Oiii] 


-0.026 


1240 ± 12 


-16.0 


1231 ± 57 


512 ± 12 


-11.4 


509 ± 26 


5007 [Oiii] 


-0.038 


3801 ± 13 


-50.1 


3762 ± 167 


1492 ± 9 


-33.0 


1479 ± 66 


5876 Hei 


-0.203 


102 ± 14 


-1.6 


97 ± 14 


96 ± 16 


-2.3 


92 ± 16 


6312 [Siii] 


-0.264 


11 ± 2 


-0.2 


10 ± 2 


16 ± 5 


-0.5 


16 ± 5 


6548 [Nil] 


-0.296 


28 ± 5 


-0.5 


26 ± 5 


61 ± 9 


-1.8 


57 ± 8 


6563 Ha 


-0.298 


3038 ± 29 


-57.0 


2797 ± 38 


3003 ± 15 


-89.3 


2807 ± 20 


6584 [Nil] 


-0.300 


136 ± 10 


-2.5 


126 ± 9 


186 ± 11 


-5.6 


174 ± 10 


6717 [Sii] 


-0.318 


277 ± 10 


-5.0 


254 ± 10 


314 ± 9 


-9.5 


292 ± 8 


6731 [Sii] 


-0.320 


188 ± 10 


-3.4 


172 ± 9 


217 ± 10 


-6.6 


202 ± 9 


7136 [Arm] 


-0.374 


112 ± 7 


-2.4 


101 ± 7 


76 ± 7 


-2.7 


70 ± 7 


7319 [On] 


-0.398 








46 ± 6 


-1.6 


42 ± 6 


7330 [On] 


-0.400 








35 ±4 


-1.2 


32 ± 4 


I(H/3)(ergseg-icm-2) 






2.70e-15 




3.19e-15 




c(H/3) 






0.12 ± 0.02 




0.10 ± O.Oi 


) 


A(V) 






0.26 ± 0.04 




0.21 ± 0.04 






NGC 6789-C 




NCC 6789-D 


A (A) 


f{A) 


F{A) 


EW{A) 


1(A) 


F(A) 


EW(A) 


1(A) 


3726 [On] 


0.322 


940 ± 13 


-29.1 


1020 ± 64 


1159 ± 26 


-15.4 


1267 ± 82 


3729 [On] 


0.322 


1377 ± 12 


-42.4 


1494 ± 92 


1712 ± 21 


-22.5 


1871 ± 116 


3835 H9 


0.299 


47 ± 9 


-1.4 


50 ± 10 








3868 [Neiii] 


0.291 


300 ± 12 


-5.2 


323 ± 23 


105 ± 9 


-0.9 


114 ± 12 


3889 Hei+H8 


0.286 


169 ± 12 


-5.5 


182 ± 17 


167 ± 20 


-2.4 


180 ± 24 


3968 [Neiii]+H7 


0.266 


74 ± 9 


-2.3 


79 ± 11 








4102 H(5 


0.229 


246 ± 16 


-5.9 


261 ± 22 


244 ± 7 


-3.2 


260 ± 16 


4340 H7 


0.157 


446 ± 7 


-12.4 


464 ± 26 


444 ± 10 


-5.3 


464 ± 27 


4363 [Oiii] 


0.149 


59 ± 3 


-1.3 


61 ± 5 


24 ± 8 


-0.2 


25 ±8 


4471 Hei 


0.115 


37 ± 5 


-0.7 


38 ± 6 


36 ± 6 


-0.3 


37 ± 7 


4686 Hen 


0.050 


26 ± 5 


-0.5 


26 ± 5 








4861 H/3 


0.000 


1000 ± 8 


-27.1 


1000 ± 47 


1000 ± 10 


-12.1 


1000 ± 47 


4959 [Oiii] 


-0.026 


1170 ± 9 


-26.9 


1163 ± 53 


607 ± 11 


-6.2 


602 ± 29 


5007 [Oiii] 


-0.038 


3471 ± 11 


-75.0 


3438 ± 153 


1750 ± 12 


-16.7 


1732 ± 78 


5876 Hei 


-0.203 


96 ± 14 


-2.8 


91 ± 14 


104 ± 20 


-1.3 


98 ± 19 


6312 [Siii] 


-0.264 


19 ± 2 


-0.6 


18 ±2 


22 ± 6 


-0.3 


20 ± 5 


6548 [Nil] 


-0.296 


51 ± 11 


-1.6 


47 ± 10 


51 ± 12 


-0.8 


48 ± 11 


6563 Ha 


-0.298 


3000 ± 23 


-92.8 


2797 ± 30 


3026 ± 16 


-46.7 


2814 ± 21 


6584 [Nil] 


-0.300 


145 ± 11 


-4.5 


135 ± 10 


149 ± 6 


-2.2 


138 ± 6 


6678 Hei 


-0.313 


27 ± 5 


-0.8 


25 ± 4 








6717 [Sii] 


-0.318 


307 ± 7 


-9.4 


285 ± 7 


300 ± 17 


-4.4 


277 ± 16 


6731 [Sii] 


-0.320 


236 ± 8 


-7.2 


219 ± 8 


206 ± 14 


-3.0 


191 ± 13 


7065 Hei 


-0.364 


23 ± 5 


-0.7 


21 ± 4 








7136 [Arm] 


-0.374 


87 ± 11 


-2.7 


80 ± 10 


80 ± 16 


-1.3 


73 ± 15 


7319 [On] 


-0.398 


44 ± 4 


-1.4 


40 ± 4 








7330 [On] 


-0.400 


42 ± 3 


-1.3 


38 ± 2 








I(H^)(ergs-icm-^) 






3.93e-15 




2.61e-15 




c(H/3) 






0.11 ± 0.02 




0.12 ± 0.05 


) 


A(V) 






0.24 ± 0.04 




0.26 ± 0.04 
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Table 2. - Continued 



A (A) 



f(A) 



NaC 6789-E 
F{A) EW(A) 1(A) 



3726 [On] 


0.322 


347 ± 51 


-30.5 


492 ± 94 


3729 [On] 


0.322 


398 ± 25 


-28.3 


563 ± 77 


3868 [Nciii] 


0.291 


62 ± 11 


-2.2 


86 ± 18 


4102 US 


0.229 


145 ± 10 


-7.6 


186 ± 25 


4340 H7 


0.157 


374 ± 9 


-11.9 


443 ± 48 


4363 [Oiii] 


0.149 


20 ± 9 


-0.6 


23 ± 11 


4471 Hei 


0.115 


36 ± 5 


-1.7 


41 ± 7 


4861 H/3 


0.000 


1000 ± 13 


-35.7 


1000 ± 93 


4959 [Oiii] 


-0.026 


842 ± 15 


-22.9 


818 ± 75 


5007 [Oiii] 


-0.038 


2358 ± 17 


-60.0 


2263 ± 201 


5876 Hci 


-0.203 


124 ± 23 


-4.3 


100 ± 20 


6548 [Nil] 


-0.296 


42 ± 14 


-1.3 


39 ± 13 


6563 Ha 


-0.298 


3090 ± 32 


-98.9 


2836 ± 42 


6584 [Nil] 


-0.300 


69 ± 12 


-2.1 


64 ± 11 


6717 [Sii] 


-0.318 


225 ± 23 


-8.4 


202 ± 21 


6731 [Sii] 


-0.320 


207 ± 18 


-8.4 


186 ± 17 


7136 [Arm] 


-0.374 


80 ± 8 


-2.8 


68 ± 6 


7319 [On] 


-0.398 


27 ± 4 


-1.2 


22 ± 3 


7330 [On] 


-0.400 


20 ± 8 


-1.0 


17 ± 6 



I(H/3)(ergs- 

c(H/3) 

A(V) 



2.94e- 
0.47 ± 
1.00 ± 



15 

0.06 

0.10 



Table 3. Electron densities and temperatures for the five knots. Densities units are in cm '^ and temperatures in 10 K. 





NGC 6789-A 


NGC 6789-B 


NGC 6789-C 


NGC 6789-D 


NGC 6789-E 


n([Sll]) 


80: 


60: 


60: 


90: 


300: 


n([Oll]) 


60: 


70: 


60: 


70: 


250: 


t([Olll]) 


1.44 ± 0.05 


1.36 ± 0.10 


1.44 ± 0.05 


1.30 ± 0.15 


1.16 ± 0.18 


t([Oii]) 


1.35 ± 0.06" 


1.18 ± 0.10 


1.41 ± 0.08 


1.21 ± 0.10° 


1.33 ± 0.25 


t([Siii])'' 


1.38 ± 0.16 


1.30 ± 0.19 


1.38 ± 0.16 


1.22 ± 0.24 


1.06 ± 0.25 



" From a relation with T([Olll]) based on photoionisation models 
^ From an empirical relation with T([Olll]) 



helium lines have negligible dependence on optical-depth ef- 
fects and the observed knots have low densities. 

The results obtained for each line and their correspond- 
ing errors are presented in Table 4. For knot C, the total 
abundance of He was calculated by adding directly the two 
ionic abundances, He/H — (He+ 4- He^"'')/H+. Also given 
in the table is the adopted value for He+/H^, a weighted 
average of the values, using the error of each line as weight. 



2.4-2 Ionic abundances from forbidden lines 

The oxygen ionic abundance ratios, 0^/H+ and O'^^/H^, 
were derived from the [On] AA 3726,3729 A and [Olll] 
AA4959, 5007 A lines, respectively using the appropriate 
electron temperature for each ion. 

The ionic abundance of nitrogen, N"'"/H"'", was derived 
from the intensities of the [Nil] AA 6548,6584 A lines assum- 
ing that T([Nll]) ^ T([Oll]). 

For sulphur, we derived S"*"/!!"^ abundances from the 
fluxes of the [Sll] emission lines at AA 6717,6731 A assuming 
that T([Sii]) « T([Oii]). Since we were not able to mea- 
sure the near infrared [Slll] AA 9069,9532 A we derived the 
S^+/H+ abundances using T([Siii]) and [Siii] A 6312 A by 
means of the equation 



12 + log{S''+/H+) = 



log 



/(6312) 



+ ■ 



I{HP) 
1.672 



-f 6.74 



t 



0.595 log t (4) 



Neon ionic abundance was estimated from the [Ne ill] 
emission line at A 3869 A. For this ion, we took the electron 
temperature of [Olll] as representative of the high excitation 
zone (Peimbert & Costero 1969). 

[Ar III] A 7136 A was the only argon line detected in the 
spectra and the abundance of Ar^"*" was calculated assuming 
that T([Arni]) « T([Sin]) (Garnett 1992). 

The ionic abundances -and their corresponding errors- 
for each observed element for the five knots are given in 
Table 5. 



3 OPTICAL AND UV PHOTOMETRY 

3.1 Hq and GALEX photometry 

We analysed the Hq image of the galaxy, retrieved from the 
Palomar/Las Campanas atlas of blue compact dwarf galax- 
ies (Gil de Paz et al. 2003). We defined elliptical apertures 
on this image for each of the five knots extracted in the spec- 
troscopic observations and we measured all the flux inside 
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Table 4. Ionic chemical abundances for helium. 



NGC 6789-A NGC 6789-B 



NGC 6789-C 



NGC 6789-D NGC 6789-E 



He+/H+ (4471) 
He+/H+ (5876) 
He+/H+ (6678) 
He+/H+ (7065) 
He+/H+ 
He2+/H+ (4686) 
He/H 



0.077 ± 0.011 
0.075 ± 0.011 



0.080 ± 0.015 
0.071 ± 0.013 



0.076 ± 0.010 0.075 ± 0.012 



0.079 ± 0.013 
0.072 ± 0.012 
0.071 ± 0.013 
0.074 ± 0.016 
0.074 ± 0.009 
0.0023 ± 0.0006 
0.076 ± 0.009 



0.076 ± 0.014 
0.076 ± 0.015 



0.082 ± 0.015 
0.073 ± 0.016 



0.076 ± 0.011 0.077 ± 0.011 



Table 5. Ionic chemical abundances derived from forbidden emission lines, ICFs" 
than helium. 



and total chemical abundances for elements heavier 



NGC 6789-A NGC 6789-B NGC 6789-C NGC 6789-D NGC 6789-E 



12 + log(0+/H+) 

12 -I- log(02+/H+) 

12 + log(0/H) 

12 -I- log(N+/H+) 

ICF(N+) 

12 + Iog(N/H) 

log(N/0) 

12 + log(S+/H+) 

12 + log(S2+/H+) 

ICF(S+ + S2+) 

12 + Iog(S/H) 

log(S/0) 

12 -I- log(Ne2+/H+) 

ICF(Ne2+) 

12 + Iog(Ne/H) 

log(Ne/0) 

12 -I- log(Ar2+/H+) 

ICF(Ar2+) 

12 + log(Ar/H) 

log(Ar/0) 



7.43 ± 0.04 
7.63 ± 0.04 
7.84 ± 0.04 
6.03 ± 0.04 

3.57 
6.58 ± 0.04 
-1.26 ± 0.06 
5.70 ± 0.03 
5.88 ± 0.16 

1.11 
6.20 ± 0.11 
-1.64 ± 0.12 
7.00 ± 0.06 

1.39 
7.14 ± 0.06 
-0.70 ± 0.07 
5.67 ± 0.08 

1.12 
5.72 ± 0.08 
-2.12 ± 0.09 



7.82 ± 0.11 
7.28 ± 0.07 
7.93 ± 0.10 
6.37 ± 0.07 

1.57 
6.57 ± 0.07 
-1.37 ± 0.12 
5.91 ± 0.07 
6.14 ± 0.23 

1.01 
6.35 ± 0.18 
-1.58 ± 0.20 
6.77 ± 0.10 

1.61 
6.98 ± 0.10 
-0.95 ± 0.14 
5.56 ± 0.12 

1.21 
5.64 ± 0.12 
-2.28 ± 0.16 



7.42 ± 0.07 
7.59 ± 0.04 
7.81 ± 0.05 

6.08 ± 0.06 
4.57 

6.74 ± 0.06 
-1.07 ± 0.08 

5.75 ± 0.05 

6.09 ± 0.16 
1.24 

6.41 ± 0.11 
-1.41 ± 0.12 
7.00 ± 0.06 

1.28 
7.11 ± 0.06 
-0.70 ± 0.08 
5.56 ± 0.09 

1.15 
5.62 ± 0.09 
-2.19 ± 0.10 



7.73 ± 0.11 
7.41 ± 0.13 
7.90 ± 0.12 
6.22 ± 0.07 

1.97 
6.51 ± 0.07 
-1.38 ± 0.14 
5.83 ± 0.07 
6.35 ± 0.30 

1.01 
6.48 ± 0.25 
-1.41 ± 0.27 
6.68 ± 0.17 

2.07 
7.00 ± 0.17 
-0.90 ± 0.21 
5.63 ± 0.17 

1.14 
5.68 ± 0.17 
-2.22 ± 0.21 



7.16 ± 0.23 
7.68 ± 0.17 
7.80 ± 0.18 
5.85 ± 0.16 

4.09 
6.50 ± 0.16 
-1.30 ± 0.24 
5.70 ± 0.14 

7.64'' 

6.58 ± 0.14 

-1.22 ± -0.25 

6.75 ± 0.20 

1.22 
6.82 ± 0.20 
-0.96 ± 0.27 
5.70 ± 0.19 

1.10 
5.74 ± 0.19 
-2.06 ± 0.22 



" ICFs estimated from tailored photoionisation models (see Section 4.2) 
>> ICF(S+) 



the elliptical apertures up to the isophote corresponding to 
the 50% of the peak of the intensity for each knot (see Fig- 
ure 1). The total fluxes and the corresponding luminosities 
at the adopted distance are listed in Table 6 with their cor- 
responding errors. The observed Ha fluxes were corrected 
in each knot for dust extinction using the values of c(H/3), 
given in Table 2. We also list in the same table the size of 
each knot, as the radius of the circular aperture whose area 
is equal to that encompassed by the elliptical aperture. As 
can be seen in Figure 1, NGC 6789 has an Hq major axis 
of 740 pc and 560 pc for its minor axis, measured using 
an elliptical aperture that encompasses the isophote at So- 
over the level of the sky background noise. The total flux 
inside this isophote (not extinction-corrected) is (2.0 ± 0.1) 



X 10 ^^ erg s ^ cm 



The Galaxy Evolution Explorer (GALEX; Martin et al. 
2005; Morrissey et al. 2005) far-ultraviolet (FUV; Kef = 
1530 A, AA = 400 A) and near-ultraviolet (NUV; Are/ = 
2310 A, AA = 1000 A) images of NGC 6789 were retrieved 
from the Nearby Galaxies Survey (NGS). We used the same 
elliptical apertures defined in the Ha image to obtain the 
FUV and NUV fluxes in each knot. In Table 6 we show for 
each knot both the FUV and NUV magnitudes once extinc- 
tion corrected together with the colour index FUV-NUV, 
all of them in AB mags. As typical sizes of the studied star- 



forming knots are lower than the GALEX point-spread func- 
tion (PSF), some aperture effects to the measured UV fluxes 
that cannot be quantified may exist. However, the character- 
ization of the knots by the corresponding colour indices are 
less affected by this effect than the total flux measurements 
in these regions, as we checked by taking different aperture 
sizes around the position of the knots. 

As can be seen, the four brightest knots (from A to 
D) have similar UV luminosities both in FUV and NUV, 
being the knot D the brightest one, although this is not 
the brightest knot in Ha. Regarding colours, knots B and C 
present redder UV colours as compared with the other three 
knots. 



3.2 WFPC2 Photometry 

Hubble Space Telescope (HST) Wide Field Planetary Cam- 
era 2 (WFPC2) images of NGC 6789 were retrieved from 
the HST archive. We discuss photometry based on the im- 
ages taken in 2000 July-September (GO-8122) through two 
continuum filters: F555W (V) and F814W (I). Each camera 
images onto a Loral 800x800 CCD which gives a plate scale 
of 0'.'046 pixel" ^ for the PC camera and O'.'lO pixel" ^ for the 
three WF cameras, with a readout noise of ~5 e" and a gain 
of 7 e^ /DN for this observations. The central star- forming 
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Table 6. Properties of the individual knots as measured in the Ho and GALEX photometry. 



NGC 6789-A NGC 6789-B NGC 6789-C NGC 6789-D NGC 6789-E 
:F(Ha) (erg/s/cm^) -13.96 ± 0.03 -13.85 ± 0.03 -13.69 ± 0.03 -14.03 ± 0.03 -14.17 ± 0.03 



log L(Ha) (erg/s) 
Radius (pc) 
FUV (mag) 
NUV (mag) 
FUV - NUV (mag) 



37.23 ± 0.02 

30 
20.00 ± 0.14 
19.86 ± 0.08 
0.14 ± 0.16 



37.34 ± 0.02 

34 
20.00 ± 0.16 
19.69 ± 0.09 
0.28 ± 0.19 



37.50 ± 0.03 

37 
20.18 ± 0.20 
19.86 ± 0.10 
0.31 ± 0.22 



37.17 ± 0.02 

36 
19.71 ± 0.15 
19.68 ± 0.09 
0.03 ± 0.18 



37.02 ± 0.01 

21 
21.24 ± 0.04 
21.18 ± 0.01 
0.06 ± 0.04 



region of NGC 6789 was centered in the PC camera in all 
images using three different pointings. The scale of the PC 
CCD at NGC 6789, for which we assumed a distance mod- 
ulus of {m~ M) = 27.80 (Drozdovsky et al. 2001), is 0.80 
pc pixel"^. Table 7 lists the details concerning the WFPC2 
data. Figure 3 shows the HST WFPC2 images of NGC 6789 
in F555W and F814W filters. The contours are drawn from 
FUV and NUV GALEX images, respectively. The aperture 
of knot E defined in the Ha image and the rectangle encom- 
passing the FOV of the PC chip of the WFPC2 data are 
also shown for reference. 

The stellar photometric analysis was performed with 
the HSTphot package (Dolphin 2000b). This package is 
specifically designed for use with HST WFPC2 images and 
uses a library of Tiny Tim (Krist 1995) undersampled PSFs 
for different locations of the star on the camera and of the 
star within the pixel, to centre the star and to find its mag- 
nitude, given in the flight system magnitude. The first step 
was to run the mask task using the HST data quality file 
(elf) to mask out the bad pixels and other image defects. 
The next step was to run crmask for cosmic ray removal. 
We used a registration factor of 0.5 and a-threshold of 3. It 
has the capability of cleaning images that are not perfectly 
aligned, and it can handle images from different filters. Af- 
ter cosmic-ray rejection, sets of images of each filter at a 
common pointing were combined into a single image, using 
the routine coadd. Since there are two images per pointing, 
we ended up with 3 images per filter. The sky computa- 
tion is made by get sky, which takes all pixels in an annulus 
around each pixel, determines the sky value and calculates 
the sky background map. These sky values are used only as 
a starting guess in the HSTphot photometry. The final step 
requires the use of the hotpixels procedure on each com- 
bined image, which uses the result from getsky and tries 
to locate and remove all hot pixels. This is an important 
step, since hot pixels can create false detections and also, 
can throw off the PSF solutions. 

The main hstphot routine was run on the images in 
the F555W and F814W bands. This task performs steUar 
PSF photometry on multiple images from difl^erent filters 
and pointings (providing the dithering pattern, it analyses 
the dithering image as a whole), including alignment and 
aperture corrections, as well as PSF modifications to cor- 
rect for errors of geometric distortion via the Holtzman et al. 
(1995) distortion correction equations and the 34th row er- 
ror, noted by Shaklan et al. (1995) (see also Anderson & 
King 1999), and correction for charge transfer inefficiency 
(Dolphin 2000a). We enabled the determination of a "local 
sky" value (option 2). 

We selected "good stars" from the hstphot output. Ob- 
ject types were classified as good star, possible unresolved 
binary, bad star, single-pixel cosmic ray or hot pixel, and ex- 
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Figure 4. CMD of NGC 6789 with average photometric uncer- 
tainties per magnitude bin. 



tended object. To ensure selecting high fidelity point sources, 
we also used the "sharpness" parameter (absolute value to 
be ^ 0.35), a measure of the quality of the fit (x^ ^ 2.5) 
and a minimum signal-to-noise ratio of 10 to reject false star 
detections in regions with structured nebulosity or artifacts. 
The final number of stars detected in the PC chip with these 
parameters in both filters were 7410. 

To quantify the completeness and systematic uncer- 
tainty of the photometry, a grid of artificial stars was gen- 
erated on a 2-dimensional CMD and distributed according 
to the flux of the images with an artificial star routine pro- 
vided by HSTphot. The parameters of the routine are the 
minimum and maximum of the measured colour and magni- 
tude. The magnitude steps used were multiple of 0.5, while 
colour steps are by 0.25. The artificial stars were distributed 
on the CMD in accordance with the number observed. Ap- 
proximately 60,000 fake stars were added in each image (at 
different trials, in order to leave the crowding conditions un- 
altered) and were given random magnitudes and colours in 
the observed range. The 50% completeness of the F555W fil- 
ter is reached at 26.9 magnitude, while for the F814W filter 
is 25.6 mag. 

The final CMD of NGC 6789 with typical photometric 
errors per magnitude bin is shown in Figure 4. 
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Table 7. Journal of HST/WFPC2 observations of NGC 6789. The images were obtained in 2000 July-September for the cycle 8 program 
GO-8122, with Rogina Schultc-Ladbeck as PI. The object of study is centered in the PC chip. 



Filter 



F555W 



F814W 



RA 


Dec 


Exposure 


(J2000) 


(J2000) 


(s) 


19:16:37.77 


+63:58:37.2 


2 X 1300 


19:16:37.73 


+63:58:37.4 


2 X 1400 


19:16:37.68 


+63:58:37.6 


2 X 1400 


19:16:37.77 


+63:58:37.2 


2 X 1300 


19:16:37.73 


+63:58:37.4 


2 X 1400 


19:16:37.68 


+63:58:37.6 


2 X 1400 



Observation ID 



U5BF0301R, U5BF0302R 
U5BF0303R, U5BF0304R 
U5BF0305R, U5BF0306R 
U5BF0401R, U5BF0402R 
U5BF0403R, U5BF0404R 
U5BF0405R, U5BF0406R 




43.0= 19'' 16'" 42.0' 41.0= 

R.A. (2000.0) 




44.0= 43.0= 19'' 16"' .12.0= 41.0= 

R.A. (2000.0) 



Figure 3. HST WFPC2 images of NGC 6789 in continuum filters F555W and F814W. The solid contours are drawn from FUV and 
NUV GALEX images, respectively. The dashed contour corresponds to the isophote at 3cr over the level of the sky background of the 
Palomar/Las Campanas Ha image. The aperture of knot E defined in the Ha image and the rectangle encompassing the FOV of the PC 
chip of the WFPC2 data are also shown for reference. North is up and east is towards the left-hand side. 



4 DISCUSSION 

4.1 Star formation history and stellar populations 

Studies to obtain information on the star formation of com- 
posite stellar systems have been proven to be very successful 
(see e.g., Aparicio et al. 1997; Hernandez et al. 2000; Dol- 
phin 2002). Any SFH recovery method relies on the assump- 
tion that a composite stellar population can be considered 
simply as the combination of SSPs, assigning a certain rel- 
ative weight to each SSP. The observed CMD is compared 
with theoretical ones created via Monte-Carlo methods for a 
variety of IMFs, binary fractions, star-formation laws, etc., 
extracting the stellar information from isochrones or stel- 
lar evolution tracks. Observed and theoretical CMDs are 
divided into boxes and converted into two-dimensional his- 
tograms of stellar density as function of colour and magni- 
tude (Hess diagrams) and compared using statistical meth- 
ods. 

For the analysis of this work, we used the StarFISH 
code'' developed by Harris & Zaritsky (2001). This code was 



successfully used in a number of cases (e.g. Harris & Zarit- 
sky 2004; Brown et al. 2006; Williams et al. 2007; Harris 
& Zaritsky 2009; Garci'a-Benito et al. 2011). Using deter- 
minations of the interstellar extinction, photometric errors, 
and distance moduli, it uses minimisation of a chi-squared- 
like statistics technique to find the linear combination of 
single-component stellar population models that best fit an 
observed CMD. The ages and metallicities of the underlying 
stellar population are characterized by the ages and metal- 
licities of the CMDs included in the best, while SFR at each 
age is provided by the weights given to the CMDs. 

The theoretical isochrones chosen for analysis were 
those of Marigo et al. (2008)^ for ages in the range 1 Myr- 
14Gyr. This data set includes models with age bins spaced 
logarithmically since the CMD changes much more rapidly 
at young ages than at old ones. We chose the two metal- 
licity values closest to the observed ones (see Section 2.4), 
namely Z = {0.001, 0.004}. Marigo et al. (2008) provide 
metallicities in the range 0.0001 ^ Z ^ 0.03. Nevertheless, 
for the sake of consistency we selected only the same values 
that were available for the STARLIGHT libraries (Bruzual 
& Chariot 2003). The photometric error and completeness 



Available at http://www.noao.edu/staff/jharris/SFH/ 
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estimates were taken directly from the results of the artifi- 
cial star experiments described in Section 3.2. We adopted 
a Salpeter IMF with a spectral index of -1.35 from 0.1-100 
Mq Salpeter (1955). This assumption is likely to be valid 
since our CMD does not contain stars with masses < 1 Mq. 
The binary fraction was set to a value of 0.5. 

Our reference value for NGC 6789 distance modulus is 
(m - M) = 27.80 ± 0.13 ± 0.18, value reported by Droz- 
dovsky et al. (2001) using the tip of the red giant branch 
distance method. The Galactic extinction is Av = 0.212, 
value provided by Schlegel et al. (1998) extinction maps for 
an area with radius equal to 5 arcmin around NGC 6789, 
remarkably close to the reddening value obtained from our 
spectroscopic data (Section 2.2). Nevertheless, we built a set 
of models to explore the space of parameters. The SFH re- 
covery is repeated for each point in the grid (m - M) vs Ay, 
and then we built a final Xmin rnap for the solutions. We 
explored the space of parameters in steps of 0.01 in distance 
modulus and extinction. To evaluate the errors of the re- 
covered solutions, we generated a series of synthetic CMDs 
using the best-fitting SFR and find a correspondence be- 
tween the Xmin and the confidence level of significance. 

Although it is tempting to perform a SFH analysis for 
each knot, the resulting individual CMDs, taken to be as the 
stars included inside the defined ellipses (see Figure 1 and 
subsection 3.1), do not contain enough stars** and, therefore, 
the large errors associated with each individual SFH prevent 
us to draw any conclusion. At any rate, the distribution of 
stars in each individual CMD is very similar for all the knots. 
Thus, we use the entire CMD of the galaxy to derive the 
global SFH. 

Figure 5 shows the overall best-fitting StarFISH solu- 
tion for the SFH of NGC 6789, located at (m - M) = 27.83 
± 0.06 and Ay = 0.64 ± 0.08, where the errors stand for 
the la confidence level. The errors in the star formation di- 
agram only allow the consideration of a few main bursts in 
the SFH. Regarding the lower metallicity (Z — 0.001), the 
solution shows a few bursty events around 1.6 Gyr and 6- 
10 Gyr, followed by a more recent burst around 100 Myr 
and a very young one during the last 12 Myr. As for the 
higher metallicity (Z — 0.004), only one event between 150- 
300 Myr is clearly seen. Although for the very young stars 
(< 12 Myr) both metallicities present a significant error, we 
take Z = 0.001 as the metallicity for the youngest event, 
as derived from the nebular analysis (see section 2.4). The 
SFR, according to StarFISH's results, was stronger around 
100 Myr, with a peak of 0.015 Mo/yr. 

In the light of these results, we used the age- metallicity 
solution derived from the CMD as a prior to constrain the 
libraries of the STARLIGHT code, used to fit the spectra 
of the five knots (see Section 2.2 for more details). In left 
panel of Figure 6 we show the age distribution of the light 
fraction obtained from STARLIGHT for each knot under 
the constrains mentioned above. Another fitting was carried 
out without any specific age-metallicity constrain in order 
to compare it with the results of the first fitting. These are 
shown in the right panel of Figure 6. As can be seen, all knots 
present a very young stellar population, with ages younger 
than 10 Myr, responsible for the ionisation of the gas, com- 
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Figure 5. Best StarFISH SFH fit derived from the HST optical 
observations of the resolved stellar populations for NGC 6789. 
The best distance modulus and extinction values are also shown. 



bined with other older populations. The presence of these 
old populations is confirmed with the detection of absorp- 
tion metal lines of Ca ll K A 3933 A and Ca ll H A 3968 A 
characteristic of old stellar populations, except for knot E, in 
which these features are marginally detected. The combina- 
tion of these two different stellar population is also obtained 
by Perez-Montero et al. (2010) in a sample of 10 BCDs using 
the same technique. In terms of mass, more than the 80% 
comes from the older population. The estimated total stellar 
mass and the fraction of the mass of the stellar population 
with an age younger than 10 Myr with respect to the total 
mass are given in Table 8. The stellar masses were corrected 
for aperture effects, with factors calculated using the ratio 
between the extinction-corrected Ha fluxes measured in the 
elliptical regions of the photometry and those measured in 
the extracted long-slit observations for each knot. In the 
same table we show the EW(H/?) values calculated once the 
contribution of the older stellar populations was removed, in 
order to use EW(H/3) to derive the properties of the ionising 
stellar population (see Section 4.2). The inner extinctions 
calculated by STARLIGHT for each fit are also listed in the 
same table. Finally, the results for the non-constrained case 
are given in the second row for each knot. The values from 
Table 8 show that those fits obtained from non-constrained 
libraries do not significantly alter some of the properties of 
the stellar populations, such as the total stellar mass and the 
correction factor for EW(H/3), even that some relevant dif- 
ferences appear in the age and metallicity distribution of the 
resulting SFHs. Therefore, we relied on the general chemi- 
cal evolution solution derived by StarFISH to constrain the 
age-metallicity libraries of STARLIGHT. At any rate, the 
general distribution of the light is similar in both cases. We 
recall that the CMD SFH (STARFISH) is derived for the 
whole galaxy, while the spectra SFH (STARLIGHT) is de- 
rived for each individual knot. 

Less significant is the difference of the fitted spectra be- 
tween both cases. For these two particular metallicities, the 
variation of the intensity of the first four Balmer lines on 
the residual spectra for each knot is well within the observa- 
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Figure 6. Histograms of tiie distribution (light fraction) of the best fit derived by STARLIGHT for the spectra of tlie five knots of NGC 
6789. The left panel corresponds to the fits using the constrains obtained by the CMD analysis, while the right panel shows the fits using 
non-constrained age-metallicity libraries. 



tional errors when using one or the other case. The change 
of age-metallicity distribution in the libraries only affects 
significantly in these two particular cases to the SFH in the 
assignments of the metallicity of the stellar populations, but 
the overall derived age distribution is very similar. 

Regarding the extinction, we find some discrepancies 
between the values found by STARLIGHT and those derived 
using Balmer decrements. However, the average extinction 
value for all knots (0.52 ± 0.16) is very close to the CMD 
best fitting value given by StarFish (0.64 ± 0.08). It should 
be noted, though, that in principle it is not expected to find 
the same values in the gas and the stellar population. 



4.2 Ionising stellar populations and 
photoionisation models 

Among the properties of the young ionising populations that 
can be derived from optical spectroscopy are their stellar 
masses and ages. The former can be derived from the Hq 
luminosity and EW(H/3) {e.g. Diaz 1998). Hq luminosity can 
also be used to estimate the SFR and the mass of ionised hy- 
drogen. The age can be estimated by using the EW(H/3) as 
compared to sequences of evolutionary synthesis models. In 
our case, the contamination of the older underlying popula- 
tion to the H/3 continuum was removed by using the analysis 
made with STARLIGHT. As can be seen in Table 8, these 
corrections are quite different between the knots, going from 
a 3% in the case of knot B to a 42% in the case of knot E. 



Table 8. Properties of the stellar populations as obtained with 
STARLIGHT for each knot of NGC 6789, including the inner 
extinction, the aperture corrected total stellar mass, the mass of 
stars younger than 10 Myr (Mio„«), and the EW(H/3) corrected 
for the contribution of the older population. The first row for 
each knot corresponds to the constrained age-metallicity libraries, 
while the second row is for the unconstrained case. 



ID 


A(V) 
(mag) 


logM. 


log Mion* 


EW(H/3), 

(A) 


Knot A 


0.4 


5.34 


3.69 


23 




0.4 


5.35 


3.81 


19 


Knot B 


0.3 


5.26 


3.16 


30 




0.7 


5.22 


3.68 


30 


Knot C 


0.7 


6.00 


4.03 


49 




0.6 


5.97 


3.95 


51 


Knot D 


0.5 


5.42 


3.79 


24 




0.5 


5.51 


3.66 


24 


Knot E 


0.7 


5.16 


2.71 


84 




0.5 


5.21 


2.47 


95 



Nevertheless, as already shown by Perez-Montero & 
Diaz (2007) and later in Perez-Montero et al. (2010) for 
other samples of BCDs, the corrected EW(H/?) are still much 
lower than the expected values for the ionisation of stellar 
clusters younger than 10 Myr. The Hi mass can be estimated 
using the CO observations of this object reported by Leroy 
et al. (2005) (marginally detected and thus, an upper limit) 
and using the Xco factor at the appropriate metallicity pro- 
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Table 9. Properties of the individual knots derived from tailored photoionisation models compared with the observed values. 





NGC 6789-A 


NGC 6789-B 


NGC 6789-C 


NGC 6789-D 


NGC 6789-E 




Mod. 


Obs. 


Mod. 


Obs. 


Mod. 


Obs. 


Mod. 


Obs. 


Mod. 


Obs. 


log L(Ha) (erg/s) 


37.32 


37.23 ± 0.02 


37.34 


37.34 ± 0.02 


37.51 


37.50 ± 0.03 


37.18 


37.17 ± 0.02 


37.03 


37.02 ± 0.01 


Radius (pc) 


7 


30 


33 


34 


22 


37 


10 


36 


5 


21 


FUV-NUV (mag) 


0.17 


0.14 ± 0.16 


0.28 


0.28 ± 0.19 


0.30 


0.31 ± 0.22 


0.22 


0.03 ± 0.18 


0.10 


0.06 ± 0.04 


I([Oii]/I(H/3) 


2.30 


2.34 ± 0.15 


3.33 


3.16 ± 0.20 


2.35 


2.51 ± 0.16 


3.05 


3.14 ± 0.20 


1.09 


1.06 ± 0.04 


I([Oin]/I{H/3) 


3.82 


3.76 ± 0.17 


1.47 


1.48 ± 0.07 


3.45 


3.44 ± 0.15 


1.72 


1.73 ± 0.08 


2.30 


2.26 ± 0.20 


12+log(0/H) 


7.93 


7.84 ± 0.04 


7.82 


7.93 ± 0.10 


7.79 


7.81 ± 0.05 


7.77 


7.90 ± 0.12 


7.64 


7.80 ± 0.18 


-EW(H/3 (A) 


21 


23 


30 


30 


48 


49 


29 


24 


80 


84 


Age (Myr) 


5.4 




4.2 




3.5 




5.2 




3.9 




Dust-to-gas ratio 


0.026 




0.34 




0.12 




0.15 




0.015 




Abs. factor (fa) 


5.34 




7.28 




8.21 




4.37 




2.23 





posed by Magrini et al. (2011) to derive the Hi column den- 
sity. For the adopted distance of NGC 6789 this gives a total 
mass of log(M_ff//M0) ~ 6. Thus, this implies that dust ab- 
sorption could be the main responsible of the disagreement 
between expected and observed EW(H/3). Thus, both L(Ha) 
and EW(H/?) must be corrected using the dust absorption 
factor (fd) before the derivation of the properties of the ion- 
ising cluster. This correction must be done even after the 
extinction correction, because it corresponds to the fraction 
of ionising photons absorbed by the dust and, therefore, not 
absorbed and re-emitted by the gas: 



Q{H) = U X Q,UH) 



(5) 



To estimate both fd and the ages of the ionising clus- 
ter for each knot, we made tailored photoionisation models 
using the code CLOUDY v.08.00 (Ferland et al. 1998). Our 
models assume a one-dimensional structure, with the central 
ionising cluster and a gas and dust geometry, and compute 
the emergent spectrum. For each knot we took as input the 
Starburst99 libraries (Leitherer et al. 1999; Vazquez & Lei- 
therer 2005) , based on stellar model atmospheres from Smith 
et al. (2002), Geneva evolutionary tracks with high stellar 
mass loss (Meynet et al. 1994), a Kroupa IMF (Kroupa 
2002) in two intervals (0.1-0.5 and 0.5-100 M©) with dif- 
ferent exponents (1.3 and 2.3 respectively), the theoretical 
wind model (Leitherer et al. 1992) and a supernova cut-ofT of 
8 Mq. We fixed the metallicity of the stellar populations to Z 
= 0.001 (« 1/12 Zq), which corresponds to the closest total 
oxygen abundance measured using the direct method in all 
the knots. We also fixed the density of particles to match the 
value derived in each knot from the [Sll] emission-line ratio 
and the chemical abundances in the gas (except oxygen) to 
match those derived using the direct method. Those species 
that were not measured were put to scale with the solar 
proportions measured by Asplund et al. (2009). Although 
not all species are expected to follow solar proportions, the 
deviation of this assumption does not affect the results of 
our models of all the other studied species. In each knot, 
we used an iterative method to find the values of the free 
parameters (oxygen abundance, dust-to-gas ratio, number 
of ionising photons, internal radius, filling factor, and age 
of the ionising cluster) that better fit the observed prop- 
erties (3727 A [Oii]/H/3, 5007 A [Oiii]/H/3, 12-hlog(0/H), 
L(Hq), and FUV-NUV colour index). The fitted properties 
along with the derived cluster age and fd in the correspond- 
ing model for each knot are shown in Table 9. The FUV- 



NUV colours were obtained by convolving the model output 
emergent spectrum for each knot with the shape of the corre- 
sponding GALEX filter, performing a direct measure of the 
resulting total flux, and finally calculating the corresponding 
colour in AB magnitudes. 

The assumptions taken for the models described above 
may be an oversimplification as they do not take into ac- 
count the three-dimensional structure of the ionized gas and 
the dust. Besides, we cannot take as input parameter the 
dust-to-gas ratio of this galaxy and possible uncertainties 
related with the adopted synthesis stellar atmospheres can 
affect the final result. 

As can be seen, the observed properties are generally 
well fitted by the models, by assuming different geometries 
and dust-to-gas ratios without noticeably changing the de- 
rived metallicity of the gas. According to our results, fd 
correlates better with the FUV-NUV index that with the 
reddening constant derived from the Balmer decrement in 
the optical spectra. For instance, knot C has the highest fd 
with the lowest c(H/3) and, in contrast, knots A and E have 
the lowest fd and the highest c(II/3). This could be indicative 
of a complex inner dust structure which is not well traced 
by the reddening constant derived from the Balmer decre- 
ment. Regarding the cluster ages, our results indicate that 
all knots have ionising cluster in a range of ages between 3.5 
Myr for knot C and 5.4 Myr for knot A. The youngest clus- 
ter appears in knot C, although knot E, which has a very low 
fd correction, has the lowest EW(H/3). The age sequence in 
the studied knots does not correspond to any specific spatial 
order. 

The total masses of the ionising clusters were also de- 
rived by using the number of ionising photons and the ages 
obtained from the respective CLOUDY models and the Star- 
burst99 libraries of the corresponding synthesis cluster at- 
mospheres. These are listed in Table 10. As can be seen, the 
stellar masses of the bursts in the four first knots are quite 
similar, while knot E has the lowest value. These masses are 
in agreement within 0.3 dex with the masses derived for the 
younger stellar population using the fitting to the optical 
spectrum made with STARLIGHT (Table 8). 

Once the fd were estimated for the five knots, it was 
possible to derive the properties that depend on the cor- 
rected total L(Hq:), including the mass of ionised hydrogen, 
the Star Formation Strengths (SFS, defined as the total mass 
of gas transformed into stars during the burst) and the SFR. 
These are listed in Table 10, along with the values for the 
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Table 10. Properties of ionising stellar populations in the knots as derived using the Q(H'') value obtained from the photoionisation 
models. 



ID 


log Q(H") 


log Mburst 


log M{Hii) 


log SFS 


SFR 




(erg/s) 


(Mo) 


(Me) 


(Mq) 


(10-3 Me/yr) 


Knot A 


49.94 


3.82 


3.45 


4.34 


1.04 


Knot B 


50.09 


3.65 


3.60 


4.32 


1.47 


Knot C 


50.31 


3.69 


3.82 


4.54 


2.44 


Knot D 


49.69 


3.52 


3.20 


4.09 


0.59 


Knot E 


49.25 


2.73 


2.76 


3.48 


0.21 


NGC 6789" 


51.49 




5.00 


5.89 


36.97 



° Values derived from the Ha image, not corrected for extinction. 



whole galaxy. The total ionised masses of the gas clouds in 
each knot were calculated using Osterbrock (1989): 



M^+=Q(ff°)-'^^ 



TleQs 



(6) 



assuming case B, n^ = 100 cm^'^, and as ~ 2.59 x 10^^'' 
cm"'^ s~^. The individual knots show masses of ionised hy- 
drogen of the order of hundreds of solar masses for knot E 
and up to ~ 6600 M© in knot C. In the case of the whole 
galaxy, the total flux described in Section 3.1 can be con- 
verted into (1.0 ± 0.1) X 10^ M0 of ionised hydrogen, com- 
parable to knot A of the Giant Extragalactic Hll Region 
NGC 5471 in MlOl (Garcia-Benito et al. 2011) or to the 
whole NGC 604 giant Hll region in M33 (Relano & Kenni- 
cutt 2009). 

It is known that the SFR determined from Ha is sensi- 
tive to several uncertainties, as those related with the extinc- 
tion and the assumed IMF. Besides, not always all the ionis- 
ing photons are absorbed, being the escape fraction of ionis- 
ing radiation from individual Hll regions (in nearby galaxies) 
between 15% and 50% (Kennicutt 1998). We applied the cal- 
ibration by Oti-Floranes & Mas- Hesse (2010)'^ which makes 
a distinction between instantaneous (IB) and extended burst 
(EB). We used our extinction and dust absorption corrected 
Hq luminosity (except for the value of the whole galaxy) 
and a Salpeter IMF in the range 0.1-100 Mq as input. For 
the IB case, the value given by the calibration is the SFS. 
The SFS values in Table 10 were calculated using the age 
provided by the calibration, namely 4, 5, and 6 Myr, closest 
to the age estimated from the photoionisation models. Ac- 
cording to this calibration, the mass of gas transformed in 
stars during the actual burst in each knot is of the order of 
tens of thousands of solar masses, except for knot E, which 
is one order of magnitude lower. 

We provide as well the SFR, for the EB case (continuous 
star formation), as a reference value. Regarding the SFR of 
the whole galaxy, the value estimated by means of the Ha 
flux is slightly higher than the peak obtained from the CMD 
(see Figure 5). The SFR for BCDs spans from a few 10~^ to 
several times 10^ Mo/yr (Hopkins et al. 2002; Zhao et al. 
2011), showing NGC 6789 not a particularly high value. 

Although our models are able to reproduce much of the 
observational information both in the optical and the UV, 
some of their limitations must be taken into account. It is 



^ Their web tool can be found at 

http://www.laeff .cab. inta-csic.es/research/sfr/ 



possible that the EW(H/3) could be modified owing to great 
fractions of leaked photons. However, comparing the mass 
of ionised gas in this galaxy [log(M/MQ) ~ 5, see Table 10] 
with the mass of neutral hydrogen does not support this sce- 
nario, but taking into account the oversimplified geometry of 
our models it cannot be ruled out completely. In this sense, 
as the long-slit does not cover entirely the gas-emission from 
the knots, the EW(H/3) are lower limits. Nonetheless, the 
aperture factors derived from the direct comparison with 
the Ha image do not lead to a satisfactory solution. 

Larger dust-to-gas ratios and fd than typical values for 
Hll regions with larger metallicity (e.g. Inoue 2001) were 
derived by the models in all knots in order to fit the ob- 
served EW(H/3) and FUV-NUV colours. However, these val- 
ues cannot be taken as representative of the whole galaxy 
but only for the analysed star-forming knots, considering 
that the fraction of ionised-to-neutral gas is quite small, as 
mentioned above, and that this galaxy does not present a 
very high IR luminosity (i.e. it does not appear in the IRAS 
catalogue). We have run other set of models with standard 
dust-to-gas ration values, but they do not fit the observables 
mentioned above. The dust-to-gas ratios could be consider- 
ably decreased by assuming different geometries of the gas 
and dust inside the Hll regions, but this approach is beyond 
the scope of this work. 

Finally, our results about the masses and ages of the 
ionising clusters must be used only as a guiding signpost, 
as the involved masses range in a regime (~ 10* Mq) 
where stochastic fluctuations in the IMF can be important 
(Cervino & Vails- Gabaud 2003; Cervifio & Luridiana 2006). 
In this sense, IMF sampling effects (Pflamm-Altenburg et al. 
2007, 2009) can also alter our results, as there can be a lower 
number of massive stars than expected for the mass range 
of these knots. 



4.3 Densities, temperatures and chemical 
abundances 

The derived densities, both for [Sll] and [On], are very sim- 
ilar for the first four knots, being lower than 100 cm~"^. In- 
terestingly, knot E, located in the outskirts of the galaxy, 
shows a density twice higher than the knots located in the 
central part of the galaxy. At any rate, these densities are 
typical of low density environments found in BCDs. 

Regarding electron temperatures, T([Oni]) was mea- 
sured in the five knots, with errors which depend on the 
quality of the spectrum in each knot going from 3% in knot 
A up to 13% in knot E. The first four star-forming regions 
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show very similar [Olll] temperatures, all within a relatively 
narrow range, between 13000 and 14300 K, while knot E 
has a lower temperature by about 2000 K. To our knowl- 
edge, there is no previous report on T([Oni]) for this galaxy 
in the literature. 

T([On]) was also measured in three of the knots (B, C, 
and E) and it was estimated in the other two knots by means 
of relations between these temperatures and the measured 
value of T([Olll]). The knot C, which is the brightest in Hq, 
has in average higher temperatures. 

We derived oxygen, nitrogen, sulphur, neon and argon 
total abundances in all the studied knots by taking the ionic 
abundances calculated in Section 2.4 and the ICFs estimated 
in our tailored photoionisation models. As in the case of tem- 
peratures, no other direct abundance determination exists in 
the literature for this galaxy. In average, NGC 6789 shows an 
oxygen abundance characteristic of the low values found in 
strong line BCDs, with values in all knots in the range of 12 
+ log(0/H) = 7.80-7.93, what is compatible with a similar 
abundance in all knots taking into account the errors. This 
behaviour, i.e. showing different star forming knots within a 
BCD close abundance values (i.e., < 0.2 dex), is also found in 
other objects using both long-slit spectroscopy (Papaderos 
et al. 2006; Cairos et al. 2009; Perez-Montero et al. 2009; 
Hagele et al. 2011) and integral field spectroscopy (Kehrig 
et al. 2008; Perez-Montero et al. 2011), or in Hll complexes 
in spiral galaxies with a similar spatial scale (Kennicutt et al. 
2003; Garcia- Benito et al. 2010). 

Since nitrogen and oxygen have different nucleosyn- 
thetic origins, their ratio is related to the chemical history of 
galaxies. Primary nitrogen synthesis is thought to occur in 
intermediate-mass stars in the CNO process during hydro- 
gen burning being hence independent of the initial heavy- 
element abundances, while secondary nitrogen production is 
expected to be produced in stars of all masses (Vila Costas 
& Edmunds 1993). At low metallicity, most part of nitrogen 
has a primary origin and a constant log(N/0) ratio is ob- 
served. However, the N/O values found in the studied knots 
of NGC 6789 range in an interval (from -1.38 to -1.26) is sen- 
sibly higher than the observed values for other BCDs or low- 
metallicity dwarf irregular galaxies (Izotov & Thuan 1999; 
Vilchez & Iglesias-Paramo 2003; van Zee & Haynes 2006; 
Perez-Montero & Contini 2009), which is around log(N/0) 
« -1.6. In the case of knot C, this ratio is even higher, with 
a value of -1.07. The cause of these relatively high N/O can 
be partially found in the corresponding ICFs derived by the 
models as the total nitrogen abundances derived from the 
assumption N/O ~ N"''/0''' give lower N/O values in knots 
A [log(N/0) = -1.41] and C [log(N/0) = -1.34], ahhough 
still higher within the errors than the typical N/O value for 
BCDs. 

On the other hand, several causes are cited in the lit- 
erature to explain this overabundance of nitrogen in low- 
metallicity environments {e.g. the pollution of the ISM by 
Wolf-Rayet stellar winds), but it is especially indicative in 
this case that we found this overabundance in the five knots, 
which present different evolutionary and excitation prop- 
erties. Therefore, the high N/O is apparently related with 
some other process affecting all the ISM of the galaxy. The 
fall of pristine gas in the galaxy could explain the N/O over- 
abundance (Koppen & Hensler 2005) and, at same time, the 



triggering of the star formation in different places of the 
galaxy with very short time intervals between them. 

The ratio of the alpha elements, sulphur, neon, and 
argon, to oxygen should be constant and show no depen- 
dence on the oxygen abundance, since all are products of 
a-processes in the same massive stars that make oxygen. 
The derived log(S/0) ratios range between -1.64 and -1.41 
in the four first knots and it is sensibly higher in knot E, 
but in this knot the derivation is much more uncertain as 
no [Slll] was measured. Although the ratio seems to increase 
from knot A to D, the average error (~ 0.18) prevents us 
from drawing any conclusion about the homogeneity. 

Regarding the logarithmic Ne/O ratio, knots A and C 
show the same value, -0.70, while for the rest of the knots 
the mean value is -0.94, with an average error of 0.07 dex 
for A and C and 0.20 dex for B, D, and E. 

The Ar/0 ratios found, in the range -2.28 to -2.06, show 
a very similar value for all knots, taking 0.17 dex as the 
average error (ranging from 0.09 up to 0.22 dex). 

Finally, the derived helium abundances are the same for 
all five knots of NGC 6789 within observational errors, and 
similar to the values found for other BCDs (Izotov & Thuan 
2004; Hagele et al. 2008). 



5 CONCLUSIONS 

In this paper, we present resolved stellar. Ha, and GALEX 
photometry, from different data archives, and WHT-ISIS op- 
tical spectroscopy of the five brightest star-forming knots of 
the nearest BCD galaxy NGC 6789 in order to study their 
SFH and metallicity. 

The spectroscopic observations of NGC 6789 were taken 
using ISIS double-arm spectrograph attached to the 4.2m 
WHT, which allowed the simultaneous analysis of the spec- 
tra, covering from 3650 up to 7500 A. From the long-slit 
spectra, we extracted and performed a detailed analysis 
of the five main star forming knots of the galaxy. Thanks 
to the measure of the [Olll] electron temperature in all 
knots and the measure of T([On]) in three of them, we pro- 
vide "direct method" abundances for oxygen, nitrogen, sul- 
phur, neon, and argon, not reported previously in the liter- 
ature. Our analysis indicates that this galaxy is metal-poor 
(12-flog(0/H) in the range 7.80 - 7.93) and chemically ho- 
mogeneous, with quite similar values of the studied species 
in all knots. At same time, all the knots present values of 
N/O which are higher than expected for the metal regime 
of this galaxy. 

We used optical data obtained from the HST archive 
to derive the SFH of NGC 6789 by means of the classical 
method of CMD reconstruction. We used the derived SFH in 
the program STARLIGHT to fit the optical spectrum and we 
compared these results with another non-constrained case. 
We corrected the emission-line measures from absorption 
and we derived the total stellar masses for each knot. Al- 
though the SFH obtained from the CMD is not foreseen by 
STARLIGHT in the non-constrained case, we checked that 
the total stellar mass and the ratio of older and younger stel- 
lar populations are not noticeably affected by the assumed 
SFH. Anyway, whenever possible, the extension of this type 
of studies for objects for which simultaneous resolved stellar 
photometry and optical spectroscopy are available will pro- 
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vide a much more accurate constrain to the age-metallicity- 
extinction degeneracy of these galaxies, as well as the com- 
parison between nebular and stellar properties. 

Finally, the ages and masses of the bursts of star for- 
mation in each knot were derived using CLOUDY tailored 
photoionisation models to fit the optical and UV photomet- 
ric properties and the observed emission-line ratios and cor- 
rected H/3 equivalent widths. We found that dust absorption 
factors correlate much better with the GALEX FUV-NUV 
colour index than with the reddening constants derived from 
the Balmer decrement, indicating a very complex inner dust 
extinction structure. The models predict for all knots instan- 
taneous bursts with ages in the range between 3 and 6 Myr. 
These ages do not follow any spatial trend in the galaxy im- 
age, so they are possibly related with the distribution of the 
galaxy in the line of sight. The dust absorption-corrected Ha 
fluxes were used to derive accurate SFRs for the individual 
knots. 

The combination of several observational and model 
techniques lead to a better and auto-consistent study of 
NGC 6789. The derivation of a non-typical metallicity evolu- 
tion using a CMD allows the use of the SFH to the subtrac- 
tion of the older stellar population using spectral fitting to 
the optical spectrum with STARLIGHT. This information, 
together with the derivation of accurate physical properties 
and ionic chemical abundances pointing to similar low 0/H 
and high N/O ratios in all knots allow the application of 
photoionisation models that predict dust-absorption factors 
fitting the GALEX colour indices and the derivation of ages 
(IBs in the range 3-6 Myr). 
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